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Abstract 

In recent years, one of the most interesting developments in quantum mechanics 
has been the construction of new exactly solvable potentials connected with the ap- 
pearance of families of exceptional orthogonal polynomials (EOP) in mathematical 
physics. In contrast with families of (Jacobi, Laguerre and Hermite) classical orthog- 
onal polynomials, which start with a constant, the EOP families begin with some 
polynomial of degree greater than or equal to one, but still form complete, orthogonal 
sets with respect to some positive-definite measure. We show how they may appear 
in the bound-state wavefunctions of some rational extensions of well-known exactly 
solvable quantum potentials. Such rational extensions are most easily constructed in 
the framework of supersymmetric quantum mechanics (SUSYQM), where they give 
rise to a new class of translationally shape invariant potentials. We review the most 
recent results in this field, which use higher-order SUSYQM. We also comment on 
some recent re-examinations of the shape invariance condition, which are independent 
of the EOP construction problem. 

1 Introduction 

As well known, the classical orthogonal polynomials (COP) of Jacobi, Laguerre and Hermite 
[1] have a lot of applications in applied mathematics and in physics. Since the early days 
of quantum mechanics, in particular, they have made their appearance in the bound-state 
wavefunctions of exactly solvable (ES) potentials. 

In this respect, the factorization method, initiated by Schrodinger [2] and developed by 
Infeld and Hull [3], has played an important role. This technique has been reformulated in 
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the framework of supersymmetric quantum mechanics (SUSYQM) [1] , also connected with 
the Darboux algorithm 

The concept of shape invariance (SI), introduced by Gendenshtein [6], has enriched the 
set of tools available, as its combination with SUSYQM has provided an integration method 
for the Schrodinger equation. In such a context, the class of translationally SI potentials 
is especially interesting because it includes all well-known ES potentials found in most 
textbooks on quantum mechanics. With the addition of some potentials depending upon 
more than two parameters [7], the classification given in [HI [H] was until recently thought 
to be complete 

In 2008, the introduction of the first families of exceptional orthogonal polynomials 
(EOP) pUl [TT] and the realization of their usefulness in constructing new ES quantum 
potentials [12], which were so far unknown translationally SI ones, produced a revival in 
this field. 

The purpose of the present contribution is to review the fast-growing branch of EOP 
and related ES potentials. 

2 New families of orthogonal polynomials and quan- 
tum mechanics 

The COP are often characterized as the polynomial solutions of a Sturm-Liouville problem 
in connection with the celebrated theorem of Bochner [T3j. According to the latter, if an 
infinite sequence of polynomials Pn{z), n = 0, 1, 2, . . . , satisfies a second-order eigenvalue 
equation 

p{z)P:{z) + q{z)P'^{z) + T{z)Pn{z) = XnPniz), n = 0, 1, 2, . . . , (1) 

then p{z), q{z) and r{z) must be polynomials of degree 2, 1 and 0, respectively. In addition, 
if the sequence is an orthogonal polynomial system, then it has to be (up to an affine 
transformation) one of the COP of Jacobi, Laguerre or Hermite. 

Many generalizations and extensions of these EOP families have been proposed, among 
which we may quote the Askey scheme of hypergeometric orthogonal polynomials and its 
g-analogue (see, e.g., [Sj), as well as algebraic deformations of SI potentials [T5| [16]. 
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More recently, the concept of an exceptional polynomial subspace and the closely re- 
lated notion of EOP were introduced [10]. As the COP, the latter are eigenfunctions of 
a second-order differential operator, but in contrast with Bochner's assumption, the first 
eigenpolynomial of the sequence is of degree m > 1, though the system of polynomials 
Pn{z), n = m, m + 1, m + 2, . . . , denoted by Xm, still forms an orthogonal and complete 
set with respect to some positive-definite measure. The Xi case was studied in detail in 
[Tot [TT]. where it was shown that there exist two such systems of Jacobi or Laguerre type. 

As an illustration, let us quote some properties of the latter |11], denoted by lI^\z), 
n = 1, 2, 3, 



-z — a — 1, 



n^\z) 



aia + 2), . . . 



z — a 



dz"^ z + a 



d 



[z + a + l)—-l 
dz 



[n 



■dz = 5n'n- 



T{n + a + 1) 



(2) 
(3) 
(4) 



and contrast them with corresponding properties [1] of conventional Laguerre polynomials 

Li")(2), n = 0, 1, 2, 



4")(z) = -z + a + l, 



d? , ^ . d' 



L^^}{z)L^::'\z)z^e-'dz = 5n', 



-n 



T{n + a + 1) 



10 n\ 
The former can actually be written as linear combinations of three of the latter. 



(5) 
(6) 
(7) 



L^:\z) = nL^:\z) - 2{n + a)L'^:l,{z) + {n + a)Li%{z). 



Both Xi families of EOP were then shown to have interesting applications in quantum 
mechanics [12]. For such a purpose, use was made of the point canonical transformation 
(PCT) method [171 [IE], which consists in mapping a Schrodinger equation containing some 
physically acceptable potential into a second-order differential equation satisfied by some 
polynomials by making changes of variable and of function. 
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In the case of equation for Laguerre-type EOP lIt'\z), the PCT method led to the 
Schrodinger equation for a rationally-extended radial oscillator (in units h = 2m = 1) 



+ VJ,ext(a;) 



i,'i\x) = E^^^\x), 0<x<oo, (9) 



where 



V^,ext(a;) = Vi{x) + Vl,rat(x), Vi{x) = -^uj^x^ + (10) 



Au 8a;(2/ + 1) 

''^'^^^"'^ " ux^ + 21 + 1 (u;x2 + 2/ + l)2- ^^^^ 
Here V/(x) is a conventional radial oscillator potential, where u and / denote the oscillator 

frequency and the angular momentum quantum number, respectively. The additional term 

Vi^rat{x) does not modify the behaviour of Vi{x) for large values of x, while for small values 

it has only an effect for / = as it changes Vo(0) = into Vo,cxt(0) = —Auj. No singularity 

being introduced on the half-line, the resulting potential Vi^ey±{,x) is a well-behaved quantum 

one. 

It turns out that V/,ext(a;) has the same bound-state spectrum as V/(x), namely 

=u;(2z/ + / + I), i/ = 0,1,2,.... (12) 
Whereas the bound-state wavefunctions of the latter, 

^^\x) oc x'+^e--^^^"L^y^\\ujx^) oc r^i{z)L^^\z), z/ = 0, 1, 2, . . . , (13) 
z = \ux'^, a = l + \, r7;(;z) = ^^^^ 
can be expressed in terms of Laguerre polynomials, those of the former, 

i^^Hx) oc oc ^n:Uz\ ^ = 0, 1, 2, . . . , (15) 

cjx^ + 2/ + 1 ^ z + a 

are given in terms of Laguerre-type Xi EOP. 

The same approach used for Jacobi-type Xi EOP led to a rational extension V^ jj ext(3^) 

of the well-known Scarf I potential Va,b(x) = [^(^4 — 1) + -B^] sec^ x — B[2A — 1) sec a; tan 
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f<x<|,0<i?<yl — llJ with similar properties [12] . 



^Note that some people prefer to use the so-cahed trigonometric Poschl- Teller (or Poschl- Teller I) po- 
tential V^ b{x) — A{A — 1) sec^ x + B{B — 1) csc^ x, Q < x < ^, related to the Scarf I by the changes of 
parameters and of variable A ^ A — B = A-\-B, a; = ^ (a; + -l). By the same transformation, VA,B,ext(a^) 
becomes an extended trigonometric Poschl- Teller potential g ^^^t {^) ■ 



Furthermore, it was proved [12] that the two new ES potentials Vi 
were translationally SI or, in other words, that their SUSY partner in a first-order SUSYQM 
approach was a potential of the same type but with different parameters, namely V;+i^ext(a;) 
and VA+i,B,ext{x), respectively (see section 5 for more details on SI). 

The discovery of this surprising property showed that SUSYQM techniques might prove 
convenient to build new examples of ES potentials related to Xi EOF, in the same spirit as 
previous works on algebraic deformations of SI potentials flbi [16] . This was confirmed by 
the construction of a rational extension VA,B,ext{x) of the generalized Foschl- Teller potential 
VaAx) = + A{A + 1)] csch^a; - B{2A + 1) cschxcothx, 0<a;<oo, fi>v4+l>l 
|19jJi In the same work, it was also pointed out that the extension procedure was not 
successful for any well-known ES potential owing to the possible appearance of a pole 
in the rational part. This was explicitly shown on the example of the Scarf II potential 
Va,b{x) = - A{A + l)]sech^x B{2A + 1) sechxtanhx, -00 < ,t < 00, A > 0. It 
was then suggested that VT symmetry might facilitate reconciling the construction of ES 
rational potentials to the requirement that these be singularity free. 

With a similar approach, the first examples of Laguerre- and Jacobi-type X2 EOF 
were built in connection with some new rationally-extended radial oscillator and Scarf I 
potentials [20j. A striking feature was that in each case there appeared two distinct X2 
EOF in contrast with the case of Xi where only one did exist. 

The question of arbitrary large m was successfully addressed by the construction of two 
distinct families of Laguerre- and Jacobi-type Xm EOF [211 123] , now labelled as type I and 
type II, and the thorough study of their properties [23i |24j. It then became clear that in 
the m = 1 special case the polynomials of the two families happened to be proportional, 
hence leaving only a single independent one. 

Later on, the X.^ EOF were shown to be obtainable through several equivalent ap- 
proaches to the SUSYQM procedure, such as the Darboux-Crum transformation 
the Darboux-Backlund one [28] and the prepotential method |29j . 



^Note that the hyperbohc Poschl- Teller (or Poschl- Teller II) potential Vj^^si^) = -A{A + l)sech^x + 
B{B — 1) csch X, < X < 00, is related to the generalized Poschl- Teller one by the changes of parameters 
and of variable A = A + B,B = B — A, x = ^x. As a consequence, Va.b, cxt{x) is connected with an 
extended hyperbolic Poschl- Teller potential § cxt{^)j which can be derived by the same transformation. 



Very recently, the Xm EOP (and related potentials) were generalized to multi-indexed 
families ^mi,m2,.--,mfe by making use of multi-step Darboux algebraic transformations [30] . 
the Crum-Adler mechanism [3T], higher-order SUSYQM |321 123] or multi-step Darboux- 
Backlund transformations 



3 Rationally-extended radial oscillator and Laguerre- 
type Xm EOP in first-order SUSYQM 

As a example, let us deal with the case of the radial oscillator potential Vi{x), defined in 
(fTOl) . in the framework of first-order SUSYQM. 

In such an approach, one considers a pair of SUSY partners \4l 

H^+) = A^A = — ^ + V^^Hx) - E, H^~^ = AA^ = — ^ + V^-\x) - E, 

dx^ dx^ ^^Q^ 

= —^ + W{x), A=-^ + W{x), V^^\x) = W^{x)tW'{x) + E, 

which intertwine with the first-order differential operators A and as AH'^^'> = H^~^A and 
^tff(-) = if^+M"!". Here W{x) is the superpotential, which can be expressed as W{x) = 
— (log0(x))' in terms of a factorization function E is the factorization energy and a 

prime denotes a derivative with respect to x. 

The factorization function is a (nodeless) seed solution of the initial Schrodinger equation 

+ V'^+\x)] (p{x) = E(j){x) (17) 



dx'^ 

with energy E smaller than or equal to the ground-state energy of V^(+). One may 
distinguish three cases. In case i, E = Eq'^\ (j){x) = tljQ^\x) and V^~^ has the same 
spectrum as except for Eq^\ which is removed. In case ii, E < (f){x) and (f)~^{x) 

are both nonnormalizable and V^~^ has exactly the same spectrum as V^'^'' (isospectral 
case). Finally, in case iii, E < eJ"^^ and 0"^ (x) is normalizable (while 4>{x) is of course 
nonnormalizable), so that V^*^") has an extra bound state at an energy E below the spectrum 
of V^W. 

Furthermore, from the intertwining relations, it results that the wavefunctions of V"*^"-* 
can be obtained from those of V^'^'^ with the same energy by applying the operator A. 



We shall be interested here in case ii (isospectral case). It turns out that for the radial 
oscillator potential Vi{x), there are two types of (nodeless) seed functions, whose energy E 
is less than -Eq^^ = uj{l + |) and whose inverse is nonnormalizable, namely 

(l)L(x) = Y](z)Li^H-z) oc x'+V^'^-"Li^''\-huJx^), (18) 



2" 



with 



x]{z) = ^3(2"+i)e^, = 2-3(2°-i)e-^, (20) 



and z, a defined in equation f |T4l) . The corresponding energies are Ej^ = —u^a + 2m + 1), 
m = 1, 2, 3, . . . , and E\l^ = —u{a — 2m — 1), m = 1, 2, . . . , / (< a), respectively. 

To get as partner potential a rationally-extended potential with a given Z, V'^~\x) = 
Vi,cxt{x) + C, where C is some additive constant, we have to start from a conventional radial 
oscillator potential V;/(x) with a different /'. The rational part of Vi^cxt{x) (see equation flTU]) ) 
can then be expressed in terms of some mth-degree polynomial in z, gm\z), as 

2" 



Vlra.t{x) = -2U { + 2Z 

9m 



(a) \ (a) 



(21) 



where a dot denotes a derivative with respect to z. According to the choice made for 0(x), 
we may distinguish the two cases 

(z) /' = /-l, = </)L,^, gt\z) = Lt-'\-z), C = -uj, 

m = 1,2,3,...; (22) 

[^^) /' = / + l, = 0j^i^^^, <?^^)(^) = lL— ^)(^), C = a;, 

m = l,2,...,/ + l(<a + l). (23) 

From the properties of Laguerre polynomials [T], it follows that V;^cxt(2;) has no pole on 
(0,oo). 

For m = 1, for instance, we get g'^\z) = z + a for type I and g'f\z) = —z — a for 
type II. Hence there is a single rationally-extended potential, as obtained in [121 120] (see 
equations ffTOl) and ffTTj) ). In contrast, for m = 2, g^\z) = + 2{a + l)z + a{a + 1)] 



or g2 '[z) = \[z + 2(a — l)z + a{a — 1)] according to the choice made. Hence there are 
two distinct extended potentials, as first observed in [20] and confirmed for all m > 2 in 
[211 [22]. 

From the wavefunctions il)[^\x) oc iriii{z)L^^ a' = /' + |, z/ = 0, 1, 2, . . . , of V^^"*"-', 
those of V^~^ with the same energy are obtained as 



ipl \x) oc A^\^\x) oc ^^^^\ y\°'\z), n = m + u, 

9 m' [z) 



(24) 



where y\^\z) is some nth-degree polynomial in z. The latter results from the application 



of a first-order differential operator Om^ on the conventional Laguerre polynomial L^f \z)^ 

y^:\z)^o^::^L^:'\z). (25) 

Here = g^? (£ - l) - ^^"^ for type I and = gl^^ {z£ + a + l)- zgt^ for type II. 
On the other hand, by directly inserting equation (p^ in the Schrodinger equation for 
(a;), one arrives at the following second-order differential equation satisfied by yi^\z), 




9 ^ ^( 

2z—-^ ] — + [z 



.(a) ..{a) 

X gm g-m 
a)^^ + z- 



9m 



(a) 
9m 



y^^\z) = {m-n)yi^\z). (26) 



The orthonormality and completeness of ipl~\x), z/ = 0, l,2,...,on the half-line imply 
that the polynomials y'^\z), = m + z/, i/ = 0, 1, 2, . . . , form an orthogonal and complete 
set with respect to the positive-definite measure z"e~^ [gm^) '^dz. According to the choice 
made for 0, such polynomials belong to the LI or L2 family of EOF and are denoted 
by L\^^{z) or L^^^^{z). In [211 1221 121], they were normalized in such a way that their 
highest- degree term is given by (— 2;)"'/[(n — m)!m!] multiplied by (—1)™ or 1, respectively. 
In conformity with the choices made in [TUl El E3], we will drop the additional factor 
(— 1)™ in case I. Then L\^^ ,^{z) and L^a,i,n{'^) both exactly reduce to l!^\z), considered in 
equations (I2|), (Ej) and (jl]). With this convention, it is easy to see that 

d 



Li 



-1) 



m—l 



iim 



dz 



- 1 



Hm 



^a,m,ni^) 



[a 



U+1 



m 



,-1 



d 



zg, 



{a) 



(27) 
(28) 
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by comparing the highest- degree terms on both sides of the equations. 

In the present approach, to get Vi^ext{x) we have had to reparametrize the conventional 
potential by considering Vi>{x) with /' 7^ /. It is worth mentioning that an isospectral 
transformation from Vi{x) to Vi^cxt{x) with the same parameter can be obtained in the 
framework of second-order SUSYQM [191 120] • 

4 Rationally-extended radial oscillator and Laguerre- 
type Xmi,m2 EOP in second-order SUSYQM 

To illustrate the construction of multi-indexed families of EOP and corresponding 
rationally-extended potentials, let us consider the case of the radial oscillator potential 
in second-order SUSYQM (SSUSY) [32]. 

In such a setting, one starts from a pair of Hamiltonians [35|, [361 EHl [391 SO] 

= -^ + V^'\x), = -^ + V^'\x), (29) 

intertwining with two second-order differential operators 

as AH^-^'^ = H^'^'^A and A!^H^'^'> = H^^^A"^. These intertwining relations imply that the 
functions p{x), q{x), V^^\x) and V^'^\x) are constrained by the relations 

2 P , I P \ 



q[x) = -p' + P - — + 



2p \2pJ 16p2 

2 2 ' 
2 , P P \ , ^ 



16p^ ' 

where c is some integration constant. Hence the knowledge of p{x) and c determines ev- 
erything else. SSUSY offers more possibilities of manipulating the initial spectrum than 
standard SUSYQM, such as deleting or creating one or two levels, moving one level or 
leaving the spectrum unchanged. Here we shall only be concerned with the last eventuality 
(isospectral transformation). 
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In the reducible case to be considered here, and A can be factorized into a product 
of two first-order differential operators, A'^ = A'^A'^ and A = AA, where A^ and A are 
expressed as in (fT6|) . while 

A^ = -^ + W{x), A=^ + W{x). (32) 
dx dx 

Then and A may be used as ladder operators in first-order SUSYQM with = A^A 
and H^-'> = AA^ as in while W{x) = -(log0(x))' and H^+'>(j){x) = 0. Simi- 

larly, from A'^ and A one may build H^'^^ = A and if^^^ = AJ^ ^ corresponding to 
some potentials 'V'^'^\x^^ factorization function ^{x) and factorization energy such that 
\/(±)(a;) = W'^[x) T W'{x) + E, W{x) = -(log0(x))' and iJ(+)0(a;) = 0. Whenever 
V^+){x) = 1/(-)(x), so that both pairs {H^+\H^-'>) and (/fW, can be glued to- 

gether, we get a reducible SSUSY system or, equivalently, a second-order parasupersym- 
metric (PSUSY) one [Ullig. 

In the isospectral case, we have to assume E < E < = and to demand 

that both (j)~^{x) and (f)~^{x) be nonnormalizable (together with (f){x) and (f){x)). 

The connection between the SSUSY and PSUSY approaches is given by the relations 

(33) 

p(x) = ^{W + W), c = E-E, 
showing that the intermediate Hamiltonian H is both partner to H^^'' and H^'^\ 

Instead of 0(x) and 0(x), we may start from two eigenfunctions 0i(x) and 02(2;) of the 
same V^^\x) with respective energies Ei and E2 (less than -Eq^^) and such that (pi^ and 
0^"*^ are nonnormalizable. Then, = 0i, = ^402 = W(0i,02)/0i, E = Ei and E' = E2, 
where W(0i,02) denotes the Wronskian of 0i(x) and 02(a;). As a consequence, 

. . ^ _ m01^ ^ (^1 - ^2)0102 . . 

^ 2>V(0l,02) 2>V(0l,02) ' ^ ^ 

which entirely determines the SSUSY partner potential through the equation 

V^^\x) = V^^\x) + Ap'{x). (35) 

In the case of the radial oscillator potential, we know that there are two types of seed 
functions f|T8l) and (fT9|) . which may be used. As the order of 0i and 02 is irrelevant as far 
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as the final potential V^'^^x) is concerned, there are three types of possibilities for the pair 
(01; 02)- In all three cases, we have to start from a potential V^'^\x) = Vii{x) with some 
different /'. The Wronskian >V(0i(x), 02(a^)) can be written in terms of some /ith-degree 
polynomial in z, g\^\z), itself expressible in terms of a Wronskian VV(/(z), (^(z)) of two 
appropriate functions of z. The results read 



-I 



[II] 



g^^\z) = w{Lt;'\-z),Lt;'\-z)) 

l^W = \4+2, 01 = 01+2,-1' <^2 = 

W{c^u<l>2)=ujx{xY+2fg^^\z), 



J — 2,7712' 



II 

l+2,m2 ' 



< mi < m2, 

mi + m2- 1; (36) 
< nil < 1712 < a + 2, 



[lit] 



W{Lt-''\z),L. 



(-a-2). 



Vl, 

W(0l,02) 



l,mi ' 



0), 

i,m2 ' 



/i = mi + 777,2 - 1; 

< mi, < m2 < a. 



(37) 



X 



gl^\z) = ^W(LL1(-^),LL7n^)) - {z + a)Ll:l{-z)L(-^^\z) 
/i = mi + m2 + 1. 



(38) 



In all three cases, provided glj^\z) does not vanish on the half-line, the SSUSY partner 
potentials can be written as 



\/(2) = Vi-2uj 



(a) 

91^ 



+ 2z 



..(a) 

91^ 

(a) 

91^ 



(a) 



-2(^1+^2) 



(39) 



where C = —2uj, 2uj or in case (i), (ii) or (iii), respectively. Hence, up to some additive 
constant, they assume essentially the same form as K'-"^^ and V'^~^ in first-order SUSYQM 
(see equation ( |2T1) ) with the only exception that gm\z) is now replaced by gl^\z). 

In the PSUSY setting, the potential in the intermediate H is some V/_i^ext, VJ+i^ext or 
V;+i,ext, as obtained in first-order SUSYQM. If we change the order of 0i and 02, then the 
intermediate potential is different although the final one remains the same. In particular, 
in case (iii), it becomes a V;_i^ext potential. 
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It is worth noting that extending the range of mi, m2 values to include mi = in case 
(i) or (ii) and mi = or m2 = in case (iii) would not lead to any new result, because 
the corresponding gli'\z) would then reduce to some gm\z) already found in first order. 
Even for the mi and m2 values listed in fl5B]) . (1571) and fl55]) . there are some possibilities of 
coincidence. As a matter of fact, the lowest-degree new gl^\z) polynomial corresponds to 
/i = 3 and is given by g^^\z) = [z^ + 3az'^ + 3{a — l){a + l)z + {a — l)a{a + 1)]/3. It may be 
obtained in case (i) or (ii) for mi = 1, m2 = 2, or else in case (iii) for mi = m2 = 1. Up to 
second order, there are therefore three distinct cubic-type extended potentials, the first two 
being those found in first order and associated with either L.^~ {—z) or 4""-') (2), the 
third one being that corresponding to the above-mentioned g^\z). Since this conclusion 
remains true when considering /cth-order SUSYQM with A; > 2, it has been conjectured 
that there are exactly /i distinct extended radial oscillator potentials (and corresponding 
EOP families) of yuth type [32J. This assertion has actually been proved for /i = 1, 2 and 3. 

Turning now ourselves to the wavefunctions ip^\x) of H'^'^\ we note that we may obtain 
them from those of H'^^\ 

(z), 1/ = 0, 1, 2, . . . , by applying the second- 
order differential operator A. They can be written as 

¥j\x)^^!^yi-\z), n = /z + z., = 0, 1, 2, . . . , (40) 

where y^\z) is some nth-degree polynomial in 2;, satisfying a second-order differential 
equation similar to equation fl26|) . but with yU and gli'\z) substituted for m and gm\z), re- 
spectively. As the polynomials considered in section 3, those obtained here form orthogonal 
and complete sets with respect to a positive-definite measure of the type z°'e~^ [gl^^) "^dz. 
We have therefore got (in general) Xmi,m2 EOP belonging to three families, which may be 
denoted by i^^;|„,,^2,n(^)' ^"mi.ma.nl^) and i^J;|i,,^2,„(2), respectively. We choose to nor- 
malize them in such a way that their highest-degree term is equal to {—z)^/[(n—^)\mi\m2^]. 

From the factorized form A = A A, it can be shown that the Xmi,m2 EOF can be 
obtained from the EOF by applying some first-order differential operator, thereby 
generalizing equations (1271) and (125]) [33j. The results read 

L'im,,mU^) = (-l)'"Hm2 -mi)-i0SX-i,ni(^)' (41) 

12 



-i,ii 



a,mi,m2,n\ ' 



where n = /i + z/, rii = mi + z/, = 0, 1, 2, . . . , and 

1 



mi ,/i 



Qrtii 



in (HTl), while 



1 



(/mi 



(42) 
(43) 

(44) 
(45) 



in (132]) or (i3]). 

The approach used in this section has been extended to A;th-order SUSYQM 
in which case one may distinguish k + 1 famihes of Xmi,m2,....mk EOP, denoted by 



-'a,mi,m2,...,Tnk,'n\ 



and associated with k + 1 possible choices FIl'^ < q < k, for the 



set of k seed functions of type I or 11. Here it is assumed that < rrii < 1712 < 



correspond to type 1, whereas < m^+i < mq^2 < ■ ■ ■ < rrik belong to type 11. The EOP 
degree n has been shown to be given by n = /i + z/, 1^ = 0, 1, 2, . . . , with 

k 

TT?,.- — 

2^'^ ' 2 



;x = - ^g(g - 1) - \(k. - q){k - q - 1) + q{k - q). (46) 

i=l 



5 Shape invariance of rationally-extended potentials 

In first-order SUSYQM, if one knows the spectrum and wavefunctions of V^^\x) (resp. 
V^-\x)), then one is able to obtain those of its partner V'^ \x) (resp. V^'^^x)) by applica- 
tion of A (resp. A^). For some special class of potentials, however, no previous knowledge 
is needed to be able to derive the spectrum and wavefunctions of both partners. This is 
the case of so-called SI potentials [6], i.e., those potentials for which the two partners are 
similar in shape and differ only in the parameters. More precisely, V^^\x] ai) are SI if 



V^-\x-ai) = V'^+\x-a2) + R{ai), 



(47) 



where oi is a set of parameters, 02 is a function of ai (i.e., 02 = /(o-i)) and the remainder 
R{ai) is independent of a; [4]. When combined with unbroken SUSYQM (case i or iii below 



13 



equation f|T7j) ). equation p7|) indeed provides an integrability condition so that the two 
potentials are solvable by purely algebraic methods. Such a condition is only a sufficient 
one, since many ES potentials that are not SI are known. 
In case i, for instance, equation (H7|) can be rewritten as 

W\x; ai) + W'{x; ai) = W^{x : as) - W'{x; as) + R{ai), (48) 

where W{x; ai) = — (log tljl^\x; ai))' in terms of the ground-state wavefunction tljl^\x; ai) 
of V^~^\x]ai). Since equation PHI) is a difference-differential equation relating the square 
of the superpotential W and its spatial derivative computed at two sets ai, as of parameter 
values, the classification of its solutions is not an easy task. 

The simplest case corresponds to translationally SI potentials for which the parameters 
ai and as are related to each other by translation (as = ai + A) [3 [HI [9] . Other cases have 
also been studied, such as SI potentials under scaling [l3l Hi] and cyclic SI potentials [15] . 
In the last cases, the potentials are only known formally as Taylor series, whereas in the 
ffist one, they are explicitly known in a closed form in terms of simple functions. 

For the conventional radial oscillator Vi{x), whose ground-state wavefunction is given by 
iPq\x) oc rii{z) with rii{z) defined in (fT4|) . the role of ai is played by / and the superpotential 
W{x;l) = Wi{x) can be written as 

Wi(x) = -ux - = -wx^. (49) 
2 X r]i 

As expected, Vi{x) - E^^^ = W^{x) - W^/(x), where = + |) = uj{a + 1). In such a 
case, equation (HHI) translates into 

W^{x) + W;{x) = Wl,,{x) - WUx) +2co= [Vi+^ix) - E^'^] + 2oj, (50) 

which illustrates the well-known SI of Vi{x), its partner being V;+i(x) (hence as = / + 1) 
and the remainder 2u being independent of x 

It turns out that the corresponding rationally-extended potentials Vi^ext{x) = Vi{x) + 
Vi,ra.t(yX), where Vi^ra.t{x) is expressed in terms of some (nonvanishing) /xth-degree polynomial 
gji'\z) as in equation f[39l) . have very similar properties 
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This is actually a consequence of the fact that their ground-state wavefunction can be 
written as iPq\x) oc rii{z)gjj^^^\z) / gjj^\z) or, in other words, that the lowest-degree EOP 
associated with g^^\z) is y^^\z) oc g']^^^\z). From iIjq \x), we indeed get the superpotential 

Wi{x) = Wi,,{x) + Wi,2{x), 

^r'^ 9l"^\ (51) 



Wuix) = Wi(x), Wioix) = -ux 



(a+l) (a) / ' 



which satisfies the equation V/^ext(a;) - = W^{x) - Wl{x), as it should be. Then, with 
ai = I, a2 = I + I and R{ai) = 2u, equation (HSjl becomes 

W^i^) + W;{x) = W^_,,{x) - Wl^,{x) +2u= [v^+i,ext(x) - EH+^'^j + 2u, (52) 

which is a counterpart of equation (13U|) . 

Since similar results have been obtained for extended potentials associated with Jacobi 
EOP, we conclude that the extensions considered here do not spoil the nice SI property of 
the conventional potentials, from which they are built. A whole new class of so far unknown 
translationally SI potentials has therefore been created. 

The appearance of such a novel family has motivated two recent re-examinations of the 
SI condition ( 148|) . independently of the EOP construction problem, in the case where the 
parameters are linked by translation. 

In the first one [IHl HT], the non-local differential equation for W has been replaced by 
two local partial differential equations. The first of them is similar to the Euler equation 
expressing momentum conservation for inviscid fluid flow in one spatial dimension, while the 
second provides a constraint helping one to determine unique solutions. In units wherein 
h is not set equal to one (in contrast with what is done here and in most other works), all 
the solutions of these equations, which do not depend explicitly on h, have been derived 
and proved to include all conventional translationally SI potentials and nothing else. The 
generalization of the algorithm to superpotentials that contain h explicitly has given rise 
to some of the extended potentials previously built in connection with EOP fl2[ [T9] . 

In the second study [48j, as in equation (ISTI) . the superpotential is separated into two 
parts, the first one being the superpotential associated with some conventional potential. 
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The second part is then divided into two pieces again and the latter are shown to satisfy 
some compatibihty condition, for which a few solutions are provided. These give rise to 
already known extended potentials |T2l [19] or slight generalizations thereof. 

Although the two last works are interesting because they have shed some light on the 
resolution of the SI condition, they have not led up to now to any really novel example 
of translationally SI rationally-extended potential (although the authors of |171 HH] claim 
the opposite). From this viewpoint, the technique based on EOF construction seems to be 
more powerful. It has the additional advantage of simultaneously providing the solutions 
of the Schrodinger equation at a low price. 

6 Final comments 

The purpose of this concluding section is to briefly review some works connected somehow 
or other with the topics of EOF and of related ES quantum potentials. 

Some applications of EOF have been considered to other problems than the resolution 
of the conventional Schrodinger equation. The inclusion of a position-dependent effective 
mass in the latter has been dealt with [IH], as well as the replacement of the real potential 
by a complex one, either PT-symmetric or not JHUSO]- The EOF have been studied in the 
framework of the quantum Hamilton- Jacobi formalism in connection with the supersym- 
metric WKB quantization condition [51j and in that of type A A/'-fold supersymmetry and 
quasi-solvability [521 [531 [M]- They have also been shown to be useful in the context of the 
Dirac equation coupled minimally or non-minimally with some external field and in that of 
the Fokker-Flanck equation [55] . 

Some generalizations of the EOF for noninteger m, which are not polynomials but 
are expressible in terms of the (confluent) hypergeometric function, have been proved to lead 
to SI potentials jSl] . Some of the latter had been studied before and shown to belong to the 
class of conditionally ES potentials, although their SI had not been recognized [SHI ISU] . 

In the course of the construction of the first examples of type I and type II Laguerre and 
Jacobi X2 EOF [20], there occurred so-called type III polynomials, which were generalized 
later on to higher m values, hence giving rise to what is sometimes referred to as the L3 
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and J3 families [281 ISll EHl EU |62] . These polynomials, however, do not qualify to be called 
EOP, since some degrees are missing from the set of them, and the corresponding extended 
potentials are not SI either. They belong instead to classes of polynomials and potentials 
often encountered (and periodically rediscovered) in applications (see, e.g., p ^ [57 | l63 | IM]). 

Finally, it is worth mentioning that families of EOP have also been constructed in con- 
nection with Wilson, Askey- Wilson, Racah and g-Racah polynomials used in the framework 
of discrete quantum mechanics [601 ESI [661 EZ] • 
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